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Abstract. We study the properties of dark matter haloes of a LCDM model in different environments. Using the distance 
of the 5th nearest neighbour as an environmental density indicator, we show that haloes in a high density environment are 
more massive, richer, have larger radii and larger velocity dispersions than haloes in a low density environment. Haloes in high 
density regions move with larger velocities, and are more spherical than haloes in low density regions. In addition, low mass 
haloes in the vicinity of the most massive haloes are themselves more massive, larger, and have larger rms velocities and larger 
3D velocities than low mass haloes far from massive haloes. The velocities of low mass haloes near massive haloes increase 
with the parent halo mass. Our results are in agreement with recent findings about environmental effects for groups and clusters 
of galaxies from deep (SDSS and LCRS) surveys. 
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1. Introduction 

Recent analyses of deep redshift surveys of galaxies, as the 
Las Campanas Redshift Survey (LCRS), the Sloan Digital Sky 
Survey (SDSS) and the two-degree-field (2dF) Galaxy Redshift 
Survey have demonstrated the presence of environmental en- 
I hancement of groups and clusters of galaxies - groups and clus- 
ters of galaxies in high density regions are richer, more mas- 
sive, more luminous and have larger velocity dispersions than 
. groups and clusters of galaxies in low density regions (Einasto 
; et al. I2()()3al l20()3bl I2()()3cl I2()03dl Ragone et al.'2004^. 

The goal of the present paper is to study the properties of 
dark matter (DM) haloes in LCDM simulations in various en- 
vironments and to compare the properties of these haloes with 
observational results for groups and clusters of galaxies. 

In the next Section we describe the numerical model and 
the halo finding procedure used in our analysis. Then we de- 
scribe the method to calculate environmental densities around 
haloes, and study the properties of haloes in different environ- 
ments. In the last two Sections we give a discussion and sum- 
mary of our results. 

The colour figures and the three-dimensional 
distribution of the simulated haloes can be 
seen at the home page of Tartu Observatory 
( |http : //www ■ aai . ee/~maret/cosmoweb . html|| , where we 
also present a catalogue of DM halos. 

2. The N-body model and the DM haloes 

For the present study we use a flat cosmological model with 
the parameters obtained from the WMAP microwave hack- 
Send offprint requests to: M. Einasto 



ground anisotropy experiment (Bennett et al. I2003I I: the total 
matter density Q.m = 0.27, the baryonic density ilb = 0.044, 
the vacuum (dark) energy density (the cosmological constant) 
Q.A = 0.73, the Hubble constant h - 0.71 (here and throughout 
this paper h is the present-day Hubble constant in units of 100 
km s"' Mpc"') and the rms mass density fluctuation parameter 
0-8 = 0.84. 

The simulations were performed using the Multi Level 
Adaptive Particle Mesh code (MLAPM, Knebe et al. 1200 U . 
This code uses adaptive mesh techniques and adaptive force 
softening; new sub-grids are created in regions where the den- 
sity exceeds a specified threshold. This scheme allows to im- 
prove considerably the spatial resolution of the particle-mesh 
code without loss in the mass resolution, similarly to the ART 
code (Kravtsov and Klypin'1997). It is known that two-particle 
relaxation may enhance the fraction of particles in clusters. 
Adaptive force softening is found to prevent mass segregation 
caused by two-body relaxation, while in other algorithms (etc. 
GADGET), based on fixed force softening, two-body relax- 
ation is more prominent (Binney and Knebe l2002t . 

We use the results of a simulation run in a cube of 
200 Mpc size, using a 256^' mesh and the same number of 
particles. Each particle has a mass of 3.57 x lO'^/i^'Mo. The 
transfer function was computed using the COSMICS code by 
E. Bertschinger (http : //arcturus .mit . edu/cosmics/i. 
For selection of dark matter haloes we used the FoF (Friends- 
of-Friends) algorithm (Zeldovich, Einasto, and Shandarin 

We shall use the linking length 0.23 in units of the mean 
particle separation, which approximately coiTesponds to the 
matter density contrast Sn/n = 80. This value was used in 
compiling the catalogue of the Las Campanas Loose Groups 
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of Galaxies, hereafter LCLG (Tucker et al. 120001 Heinamaki 
et al. I2003> . It is substantially lower than the spherical col- 
lapse model prediction for virialized objects, Snin - 178n,„"^ 
(White, Efstathiou and Frenk 1993 1. 

Too large a linking length may lead to overestimation of 
the masses of the simulated groups since particles outside of 
the virialized core are included in the groups. To avoid pos- 
sible effects of unbound groups in simulations, we used the 
virial condition - Ekml\Epoi\ < 0.5 {Ep,,! is the potential 
energy and Ekin the kinetic energy of a group) for groups to be 
included in our final group catalogue. Moreover, to avoid too 
small groups we choose only groups which included more than 
100 dark matter particles, imposing a minimum halo mass of 
3.57 X IO'^/j^'Mq. The catalogue contained 5355 dai'k matter 
haloes, and after applying the virial condition, 5219 dark matter 
haloes remained in our final catalogue. 

Fig-Eshows the cumulative mass function for dark matter 
haloes, compared with observational data (LCLG). The mass 
function of the LCLG (with Poisson error bars) is shown by 
filled circles. For a more detailed description of selection ef- 
fects in the LCLG sample and of the determination of the 
LCLG mass function see Tucker et al. (2000 1 and Heinamaki 
et al. J2003> . The solid line shows the mass function of our 
simulations, with halo masses determined by the sum of the 
particle masses in the halo. We see that the mass functions of 
the LCLGs and of the simulation are rather similar. Less mas- 
sive haloes in our simulation correspond to galaxy groups; the 
most massive haloes correspond to rich clusters of galaxies. 




Fig. 1. The observed mass function and the simulated mass 
function for haloes; the filled circles show the LCLG result and 
the solid line describes our simulations. 



3. Haloes in different environments 

3.1. Environmental densities around haloes 

In the present work we determine spatial densities around 
haloes, using the distance to the halo's 5th neighbour halo. We 
shall use this distance as an environmental parameter. When 
we use this distance as the density indicator, we always see the 
characteristic scales of the features we study. This method has 



been used before to find environmental (projected or spatial) 
densities around galaxies (Caon and Einasto 1995, Goto et al. 
2003 1, around clusters (Schuecker et al. i2001i Einasto et al. 
2004a I and even around superclusters of galaxies (Einasto et 
al. [T997t . 

Using the 5th neighbour distance we shall divide the haloes 
into four environment regions as described in Table ^ The 
neighbour distance intervals are chosen to have approximately 
the same number of haloes in each region. The region H23.1 
has the closest neighbours and the highest environmental den- 
sities, the class H23.4 has the most distant neighbours and the 
lowest environmental densities. 

3.2. Two richness classes of haloes 

In Einasto et al. J2003c> we showed that the properties of the 
LCLGs in the neighbourhood of rich clusters are enhanced 
in comparison with the properties of the LCLGs in average. 
Moreover, in Einasto et al. ( 2003d i we showed that the LCLGs 
in superclusters are more massive than loose groups of galaxies 
that do not belong to superclusters. 

To test this result using LCDM haloes, we divided our sam- 
ple of LCDM haloes into two richness classes: the most mas- 
sive haloes that correspond to rich clusters of galaxies, and 
less massive haloes that correspond to loose groups of galaxies. 
Here we used the mass threshold M > 6.3 x lO'^/z 'M© (see 
Heinamaki et al. 12003). Einasto et al. ( I2003ct also used this 
threshold to determine the population of Abell-class clusters 
among the LCLGs. 

Once two subsamples with different masses were deter- 
mined, we found for the most massive haloes their 5 clos- 
est neighbours among less massive haloes (with masses M < 
6.3 X 10'^/;"'Mq). Then we divided the population of low 
mass neighbours into the populations of close neighbours (CN) 
and distant neighbours (DN), using a threshold distance of 6 
/j ' Mpc. This distance was used in Einasto et al. ( I2003ct to 
define the populations of loose groups around rich clusters of 
galaxies. 

4. The properties of haloes in different 
environment 

4.1. Haloes in high and low density regions 

The spatial distribution of haloes in our L - 200 /i"' Mpc sim- 
ulation box is shown in Fig.|2 (a 50 /z"' Mpc thick slice in the 
Z-direction; color figures of all slices are given in our web page 
together with animations). This figure shows haloes located in 
regions of different densities, as quantified by the 5th nearest 
neighbour distance. In general, in high density regions haloes 
are bigger and more massive than haloes located in low den- 
sity regions. The velocities of haloes in high density regions 
are also larger than the velocities of haloes in low density re- 
gions. The velocity vectors of haloes are directed toward rich 
filaments and massive haloes - mass flows from voids into fil- 
aments, which act as attractors. In extreme void regions, where 
the haloes are most isolated, the haloes almost do not move; 
they also have low masses. 
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Fig. 2. The spatial distribution of haloes in regions of different density in a slice of thickness of 50 /; ' Mpc (in the Z-direction); the 
Z interval is 150-200 /z ' Mpc. The left panel shows the spatial distribution of haloes, the right panel shows the velocity vectors 
of haloes. Filled large circles denote haloes in the highest density environment (sample H23.1 in Table filled small diamonds 
stand for the sample H23.2, open circles for the sample H23.3, and crosses denote haloes in the lowest density environment, 
sample H23.4. 

Table 1. Median and upper quartile (in parentheses) values of halo properties in high and low density environments. 



Sample 


A'halo 




£»2(A'5) 


A'p 


log My„ 


R. 


cr 


V 


E 


Nm 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


H23.1 


1315 


0.0 


7.9 


409 (992) 


12.70(13.17) 


0.40 (0.51) 


240 (380) 


480 (640) 


0.50 (0.63) 


101 


H23.2 


1287 


7.9 


9.7 


357 (816) 


12.59 (13.04) 


0.38 (0.50) 


210 (340) 


445 (585) 


0.51 (0.65) 


66 


H23.3 


1340 


9.7 


12.2 


309 (651) 


12.49 (12.92) 


0.37 (0.46) 


190 (300) 


400 (535) 


0.50 (0.64) 


37 


H23.4 


1280 


12.2 




247 (463) 


12.32(12.70) 


0.36 (0.44) 


160 (250) 


345 (470) 


0.52 (0.66) 


18 



The columns in the Table are as follows: 

column 1 : the halo sample number, 

column 2: the number of haloes in the sample, 

column 3: the smallest 5th nearest neighbour distance, 

column 4: the largest 5th nearest neighbour distance, 

column 5: the number of particles in a halo, 

column 6: the virial mass of a halo (log), M,,> = cr-i?,.„ /G, 

column 7: the virial radius of a halo (for equal mass particles), i?,-i> = ^"/(2oa/>i ~)' 

' ' 'ij 

column 8: the rms velocity of halo particles, cr = -^2; "'f IN, where the squared velocity deviation = (v,, - v„.)^ + (v,_j, - v^^,)^ + (v« - v„)^ is 
relative to the mean velocity of the halo and v„, v^., Vcz are the components of the mean velocity, 
column 9: the velocity of a halo, vi,aio = yj'^'cx + ^ly + 

column 10: the eccentricity of a halo, £ = 1 - cc/aa, where cc is the semi-minor axis and aa is the semi-major axis, 
column 11: the number of haloes with masses M > 6.3 x IO'^/j^'Mq. 



Let us now analyse the properties of haloes in various en- 
vironments in more detail. We made a series of scatter plots, in 
which we plotted the number of particles in haloes, halo's virial 
radii, masses and velocity dispersions against the distance of 
the 5th nearest neighbour of the halo, and calculated cumula- 
tive distributions of these parameters (the number of particles in 
haloes, halo's virial radii, masses and velocity dispersions) for 
haloes from different environments, and the median and upper 
quartile values of the properties of haloes from high and low 



density regions. We have also estimated the statistical signifi- 
cance of our results, using the Kolmogorov-Smirnov test. 

We compare haloes from different environments in Tabled 
The columns show median and upper quartile (in parentheses) 
values of halo properties in high and low density environments. 
In order to not to overcrowd the paper with figures, we show 
only one scatter plot, for the halo richness (Fig.|3}- 

Fig. 13 and Table ^ show that the properties of haloes in 
high density regions are environmentally enhanced: they are 
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Fig. 4. The cumulative distributions of various parameters of haloes in different environments. Upper left panel: cumulative 
richness. Upper right panel: the mass functions. Lower left panel: the cumulative virial radius. Lower right panel: the cumulative 
velocity. Solid line - the supercluster region H23.1, dotted and dashed lines - the filament regions H23.2 and H23.3, and dot- 
dashed line - the void region H23.4. 
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Fig. 3. The number of particles in haloes versus the distance of 
the 5th nearest neighbour of the halo. 

richer, larger, more massive, and have larger velocity disper- 
sions than haloes in low density regions. The haloes in high 
density regions have larger velocities than haloes in a low den- 
sity environment. In the lowest density regions (low density 
filaments and voids) haloes almost do not move. This is seen 
also in Fig.|2l 



The cumulative mass function (MF) is defined as the num- 
ber density of clusters/groups above a given mass M, n(> M). 
We used the same definition also for other properties of the 
haloes. Thus we get the cumulative richness distribution, the 
cumulative virial radius distribution, the cumulative rms veloc- 
ity distribution, and the cumulative peculiar velocity distribu- 
tion. We plot these functions for haloes in different environ- 
ments in Fig. 12 These figures and Table[T]show that haloes in 
a high density environment are richer, they have larger masses, 
larger radii, larger rms velocities and and larger peculiar veloc- 
ities than haloes in a low density environment. For example, at 
the density level 10^^/j^Mpc"^ the H23.1 clusters are about 3 
times richer than the H23.4 clusters. The cluster class haloes 
(with mass M > 6.3 x IO'-^/z^'Mq) are preferentially located in 
high density regions. A few cluster class haloes, however, are 
located in lower density environments. This is due do the selec- 
tion procedure of the haloes - the most massive haloes do not 
have very close neighbours (they are halo's members), thus the 
local density is suppressed in the neighbourhood of the halo. 
This is seen also in colour figures in our web page. 

We have determined also the eccentricities of DM-haloes. 
Fig. |5] and Table [2 hint at the trend that the eccentricities of 
haloes in high density regions seem to be smaller than the ec- 
centricities of haloes in low density regions. Haloes in high 
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Fig. 5. Left panel: the cumulative distributions of eccentricities for simulated haloes in different environments. Solid line — 
H23.1, dotted line — H23.2, dashed line — H23.3, and dot - dashed line — H23.4. Right panel: The eccentricity of haloes versus 
the richness of haloes, for all environments. 



density regions are more spherical than haloes in low density 
regions, an evidence of a higher degree of virilization. 

Fig- 13 shows also, that the differences between the eccen- 
tricities of haloes from different environments are smaller than 
the differences between other parameters (distributions on the 
Fig. |5] left panel, for samples H23.2 and H23.3 even practi- 
cally overlap). The reason for that can be seen in the right panel 
of Fig.|5l haloes of medium richness (about 200 - 1000 parti- 
cles, mass < 3.57 x IO'-'Mq) have smaller eccentricities. This 
is probably due to several factors: haloes of medium richness 
are more evolved than very poor haloes, this makes them more 
spherical. The shape of the richest haloes is more elongated due 
to the influence of the large scale structure around them (fila- 
ments). The environmental dependence of the shapes of haloes 
is affected by these factors and is thus weaker than other envi- 
ronmental dependences. 

The Kolmogorov-Smirnov test shows that the differences 
between the parameters of haloes from high and low density 
regions are statistically significant at the 99% significance level 
for all parameters except eccentricities. The differences be- 
tween eccentricities of haloes in the highest and the lowest den- 
sity regions are statistically significant at the 75% significance 
level according to the Kolmogorov-Smirnov test. 

4.2. Properties of haloes in tine neiglibourliood of tine 
most massive haloes 

In Table|2lwe give the median and upper quartile values of the 
parameters of the halo populations CN and DN (close neigh- 
bours and distant neighbours, see Sect. 13. 2> . In Fig.|6lwe show 
the cumulative richness distributions, the mass functions, the 
cumulative virial radius distributions, and the cumulative pe- 
culiar velocity distributions for low mass haloes around high 
mass haloes. In Table |2 we also give the statistical signifi- 
cances of the differences between the samples, according to the 
Kolmogorov-Smirnov test. 



Fig. |6l and Table |2] show that the low mass haloes in the 
close neighbourhood of high mass haloes are themselves richer, 
more massive, have larger rms velocities and move faster than 
the low mass haloes far from massive haloes. The eccentricities 
of haloes in the neighbourhood of massive haloes are smaller 
than the eccentricities of haloes farther away from massive 
haloes. 

We also checked whether the properties of low mass haloes 
in the vicinity of massive haloes depend on the masses of high 
mass haloes. Our calculations show that while the intrinsic 
properties of low mass haloes do not seem to depend on the 
mass of the parent massive halo, peculiar velocities of low mass 
haloes in the vicinity of massive halo increase with the parent 
halo mass (Fig.0 left panel) - an indication of the gravitational 
influence of massive haloes. 

We obtain similar results if we use summed particle masses 
instead of virial masses for the massive haloes. 

Our calculations show that if we find all low mass neigh- 
bours around high mass haloes within a 6 /z ' Mpc sphere, (not 
just 5 closest neighbours), then we obtain almost the same pop- 
ulation (with about 470 member haloes) of low mass neigh- 
bours. Thus we used the five closest neighbours, in accordance 
with our definition of environmental densities around haloes. 

5. Discussion 

5.1. Comparison with observations 

Already early studies of galaxies in different environments 
showed the presence of strong morphology-density and 
luminosity-density relations (Dressier 1 980k Postman & Geller 
1984 Einasto & Einasto 1987 Einasto fTWTal and II 99 1 bl Mo 
et al. 1992 1. Presently, these relations are studied intensively, 
using new deep surveys of galaxies (Norberg et al. i2001»i2002l 
Croton et al. 2004 Kauffmann et aL EUO?! Balogh et al. lMMI 
and Blanton et al. 2004 among others). 

However, the study of the properties of galaxy groups 
and clusters in various environments is only in the beginning. 
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Fig. 6. The cumulative distributions of various parameters for low mass haloes in the vicinity of the most massive haloes. Upper 
left panel: cumulative richness. Upper right panel: the mass functions. Lower left panel: the cumulative rms velocity. Lower right 
panel: the cumulative peculiar velocity. Solid line — close neighbours, the sample CN, dashed line — distant neighbours, the 
sample DN. 



Table 2. Median and upper quartile (in parentheses) values of the properties of low mass haloes in the neighbourhood of massive 
haloes . 



Sample 




Di(N5) 


£»2(A'5) 




log Mvir 




O" 


V 


E 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


CN 


450 


0.0 


6.0 


395 (774) 


12.73 (13.07) 


0.40 (0.50) 


240 (340) 


600 (770) 


0.51 (0.62) 


DN 


368 


6.0 




339 (711) 


12.61 (13.00) 


0.40 (0.50) 


220 (320) 


550 (690) 


0.50 (0.65) 


KS 








80% 


98% 


85% 


97% 


99% 


70% 



The columns are as in TableQ The last line shows the significance levels of the difference of the distributions, by the Kolmogorov-Smimov 
test. 



Observations show that groups and clusters of galaxies in high 
density regions, in the vicinity of rich clusters of galaxies and in 
superclusters, are more massive and luminous and have larger 
velocity dispersions than loose groups in average (Einasto et 
al . '2003 a" ^O Olbl I2003cl l2003dl K004al 2004b). This enhance- 
ment extends to scales up to about 15-20 /z ' Mpc around rich 
clusters. 

Ragone et al. (2004) determined groups of galaxies for the 
2dF Galaxy Redshift Survey and studied the properties of these 
groups in the vicinity of rich clusters. This study confirms the 
results obtained in Einasto et al. ( 2003c i: groups in the vicinity 
of rich clusters of galaxies are themselves also richer and more 
massive than groups in average. 



In Einasto et al. ( 2003d i we studied the environmental en- 
hancement of the LCLGs in superclusters. In particular, we 
calculated for the LCLGs in superclusters the distance to the 
nearest Abell cluster in supercluster and studied the properties 
of the LCLGs as a function of the distance to the nearest rich 
cluster in supercluster. In analogy to that test we plot the dis- 
tribution of rms velocities of low mass haloes in the vicinity of 
the massive haloes (Fig. El right panel). Note the similarity of 
Fig. El right panel, and Fig. 5 in Einasto et al. (2003d): groups 
close to rich clusters have larger rms velocities than groups far 
from rich clusters in a supercluster. 

Several studies of the correlation function of nearby groups 
of galaxies show that groups of higher mass are more strongly 
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Fig. 7. Left: peculiar velocities (in km/s)nof low mass haloes around high mass haloes versus the mass of the high mass halo. 
Right: the rms velocities (in km/s) of low mass haloes around high mass haloes versus the distance of the low mass halo from the 
high mass halo. 



clustered than groups on average (Giuricin et al. 120011 Girardi 
et al. 120001 and Merchan et al.|2000l. Stronger clustering is an 
indication that these groups could be located in high density 
regions (Einasto et al. .l997j . 

Plionis, Basilakos and Tovmassian ('2004'^ showed that the 
shapes of poor groups of galaxies depend on their richness, 
poorer groups being more elongated than richer groups and 
clusters of galaxies. Observations show that very rich clus- 
ters of galaxies also are elongated (see references in Plionis, 
Basilakos and Tovmassian I2004> . This is in accordance with 
our results, showing that the poorest haloes have larger ec- 
centricities than haloes of medium richness, and the most 
rich haloes again have larger eccentricities than medium rich 
haloes. 

Several studies of clusters of galaxies have provided ev- 
idence that properties of rich clusters depend on their large 
scale environment (Novikov et al. 119991 Chambers et al. 120021 
Einasto et al. 2001 Schuecker et al. 120011 and Plionis and 
Basilakos .2002j up to a distance of about 20 /i"' Mpc. This 
distance is close to the so-called "pancake scale" (Melott & 
Shandarin 1993 1, and corresponds to the mean thickness of su- 
perclusters (Einasto et al. . 1994.. 1997. and Jaaniste et al. .l998) . 

Our results about the properties of haloes in different en- 
vironments are in accordance with these observational results. 
The LCDM model haloes in high density regions are richer, 
more massive etc than haloes in low density environment. In 
particular, low mass haloes around the most massive haloes are 
themselves also more massive and have larger velocities than 
haloes farther away from massive haloes. 

Our study shows that haloes move toward massive systems. 
In the vicinity of massive systems the velocities of haloes are 
larger than in low density regions (voids). This is in accordance 
with observational studies of velocity fields of galaxy clusters 
which show that clusters are moving toward superclusters of 
galaxies (Hudson et al. .2004. and references therein). 



5.2. Comparison with simulations 

Ragone et al. J2004t determined groups of galaxies from simu- 
lations by the Virgo Consortium, and studied the properties of 
these groups in the vicinity of rich clusters. This study con- 
firmed the results obtained using observational data: groups 
in the vicinity of rich clusters of galaxies are themselves also 
richer and more massive than groups in average. 

Ragone et al. ('2004'( also found that the enhancement of 
the properties of low mass haloes in the vicinity of massive 
haloes is a strong function of the mass of these massive haloes, 
while we did not detect such a strong dependence (except for 
3D velocities). This difference is probably due to the fact that in 
this test, Ragone et al. (|2P04 1 used different host samples with 
increasing mass limit for host haloes, while we did not change 
the host sample. 

Suhhonenko J2003> has demonstrated using diff^erent N- 
body simulations that in simulated superclusters more massive 
clusters are located in the central regions of superclusters. 

Gotdober et al. (l2002l l and Faltenbacher et al. ( l2002l 
analysed high-resolution simulations of formation of galax- 
ies, groups, and clusters, and found a significant enhancement 
of the mass of haloes in the environment of other haloes. 
This effect is especially significant at scales below 10 /i^'Mpc. 
Therefore, environmental enhancement of the halo mass is a 
direct evidence for the process of the hierarchical formation 
of galaxy and cluster haloes. Halos inside clusters formed ear- 
lier than haloes in low density filaments and are more evolved 
(Gottlober, Klypin and Kravtsov lMJTI Gottlober et al UM^ . 

Kasun and Evrard ("2004*1 found that the shapes of LCDM 
haloes do not depend on the large scale (supercluster) environ- 
ment of haloes. However, the small hint of the such correlation 
seen in our study may ensue from the different definitions of 
the environmental density. 

The finding that the richest haloes have larger eccentricities 
is in accordance with the recent results by Hopkins, Bahcall and 
Bode ( 2004i) who also found that clusters are aligned with one 
another, in accordance with several observational studies (see 
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references in Hopkins, Bahcall and Bode 2004i and the results 
by Colberg et al. ( 1999 1 who showed that the formation of a 
cluster in simulations is governed by its surrounding large scale 
structure, and the properties of clusters are linked to the prop- 
erties of structures around them. More massive clusters have 
more filaments than less massive clusters (Colberg, Krughoff 
and Connolly I2004at . This finding is in agreement with our 
result about the environmental dependencies of halo's proper- 
ties. Moreover, observations show that rich superclusters usu- 
ally have a multi-branch structure and these superclusters con- 
tain more luminous clusters than poor superclusters (Einasto et 
al.|2003b*. 

Bahcall, Gramann and Cen and Colberg et al. ( EUnOt 

found evidence that rich clusters have larger peculiar velocities 
than poor clusters, and clusters in high density regions have 
larger peculiar velocities than clusters in low density regions. 

In supercluster-void network haloes move toward density 
enhancements, where the velocities of haloes are larger than 
velocities of haloes in a low density environment: rich systems 
become richer and the fraction of matter in voids decreases 
(see also Einasto et al. '1994'!. Halos in the lowest density re- 
gions have the smallest masses in agreement with the recent 
results about the mass function of haloes in voids (Colberg et 
al.|2004bV 

Our present paper extends the earlier findings about the en- 
vironmental enhancement of haloes and shows that haloes in 
higher density environment, in general, are richer, more mas- 
sive and have larger velocity dispersions than haloes in low 
density environment. 

A possible explanation of the environmental enhancement 
of galaxy systems in high-density regions was investigated by 
Einasto et al. (I2004al and I2004bl . These studies show that dy- 
namical evolution in high-density regions is determined by the 
high overall mean density that speeds up the clustering of par- 
ticles. Therefore, in high density regions clustering starts early 
and continues until the present. The haloes that populate high 
density regions are themselves also richer, more massive and 
have larger velocities than the haloes in low density regions. In 
low density filaments that cross voids, as well as in the outer 
low density regions of high density systems, the mean density 
decreases and thus the evolution is slow, and in these regions 
haloes themselves are also poor, less massive and have small 
velocities. 

6. Conclusions 

We studied the properties of dark matter haloes in different en- 
vironments. Our main results are as follows. 

- We composed a catalogue of haloes found using the 
MLAPM code and the FoF algorithm with a linking length 
of 0.23 in the units of mean particle separation. This cata- 
logue is available from our web pages. 

- Our analysis shows that haloes in a high density environ- 
ment are richer, more massive, have larger virial radii, have 
larger rms velocities, have larger peculiar velocities and are 
more spherical than haloes in a low density environment. 



- Low mass haloes in the vicinity of high mass haloes are 
themselves richer, more massive, slightly more spherical 
and have larger rms velocity dispersions than low mass 
haloes farther away from high mass haloes. The larger the 
mass of parent halo, the larger are the velocities of low mass 
haloes in the vicinity of high mass haloes. 

- Our study indicates the importance of the role of a high den- 
sity environment, which affects the properties (formation 
and evolution) of galaxy systems. In high density regions 
haloes formed earlier, and are more evolved than haloes in 
low density regions, in accordance with the scenario of the 
hierarchical formation of the structure in the Universe. 
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